Time Resolved Correlation measurements of temporally heterogeneous dynamics 
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Time Resolved Correlation (TRC) is a recently introduced light scattering technique that allows to 
detect and quantify dynamic heterogeneities. The technique is based on the analysis of the temporal 
evolution of the speckle pattern generated by the light scattered by a sample, which is quantified by 
ci{t, t), the degree of correlation between speckle images recorded at time t and t + r. Heterogeneous 
dynamics results in significant fluctuations of ci{t,T) with time t. We describe how to optimize 
TRC measurements and how to detect and avoid possible artifacts. The statistical properties of the 
fluctuations of c/ are analyzed by studying their variance, probability distribution function, and time 
autocorrelation function. We show that these quantities are affected by a noise contribution due to 
the finite number A'' of detected speckles. We propose and demonstrate a method to correct for the 
noise contribution, based on a A'' ^ oo extrapolation scheme. Examples from both homogeneous 
and heterogeneous dynamics are provided. Connections with recent numerical and analytical works 
on heterogeneous glassy dynamics are briefly discussed. 

PACS numbers: 82.70.-y, 42.30.Ms, 64.70.Pf, 05.40.-a 



I. INTRODUCTION 

Soft glassy systems such as concentrated colloidal sus- 
pensions, emulsions, surfactant phases, gels and foams 
exhibit very slow and unusual dynamics, characterized by 
non-exponential relaxations of correlation and response 
functions, which often depend on sample history and may 
be heterogeneous both in space and time . These phe- 
nomena have attracted a great interest, in part due to 
their "universal" character. Examples of unifying de- 
scriptions are the mode coupling theory of the stationary 
average dynamics of concentrated suspensions ofparti- 
cles with both repulsive or attractive interactions 
or, at a more qualitative level, the concept of jamming, 
which rationalizes the fluid-to-solid transition in a wide 
range of systems and experimental configurations l^. 
Additionally, soft glassy materials often exhibit intrigu- 
ing similarities with hard condensed matter glasses, such 
as the dependence of the dynamics on sample history 
(agingphenomena |^ Q, rejuvenation and memory ef- 
fects Ja H) or the presence of dynamical heterogeneity 

(Ca 111 111). 

Light scattering is a popular means to measure the av- 
erage dynamics. In a traditional dynamic light scatter- 
ing experiment, one measures the normalized time auto- 
correlation function of the scattered intensity: (72 (t) = 

2 

Ii{t)Ii{t + T)/Ii{t) , where the average • • • is over time, 

t, and Ii{t) is the intensity collected by a single detector. 
The intensity autocorrelation function provides quanti- 
tative information on the dynamics of the sample; the 
way this information is extracted depends on the exper- 
imental configuration. In single scattering experiments. 



g2{T) is related to the intermediate scattering function, 
/(r), via the Siegert relation: /(r) 
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where (3 is the coherence factor that depends on the size 
ratio between the speckle — or coherence area — and the 
detector In the opposite hmit of strong multiple 

scattering, the diffusing-wave spectroscopy (DWS) for- 
malism jl5l | allows the particle mean square displacement 
to be calculated from 172 (t) — 1, provided that the dynam- 
ics be spatially and temporally homogeneous. For glassy 
systems the average over time required to compute 32 (t) 
may become in practice unfeasible, either because the 
dynamics is so slow that prohibitively long experiments 
would be required to accumulate a satisfactory statistics, 
or because the dynamics may be non-stationary, e.g. for 
aging systems. Various schemes have been introduced 
to address this issue, most of them based on the idea of 
measuring (72 (t) for many independent speckles and av- 
eraging the intensity correlation function not only over 
time, but also over distinct speckles. This can be done 
either sequentially -e.g. by slowly rotating the sample so 
as to illuminate a single detector with different speckles 
(interleave [l6j or echo 17] methods)- or in parallel, e.g. 
by using the pixels of a charge-coupled device camera 
(CC D) as independent detectors (multispeckle method 

111 El 13). 

These techniques drastically reduce the required time 
average and thus extend the applicability of light scatter- 
ing to glassy systems. Similarly to traditional light scat- 
tering measurements, however, they provide information 
only on the average dynamics, not on its fluctuations. 
However, recent theoretical and simulation works sug- 
gest that dynamic heterogeneity is a key feature of the 
slow dynamics in glassy systems. In view of the very 
limited number of experiments that directly test this be- 
havior on soft glasses [lol ITU , new experimental tools 



2 



that access fluctuations of the dynamics are needed. We 
have recently introduced the Time Resolved Correlation 
(TRC) scheme jl^, a method that allows temporally 
heterogeneous dynamics to be investigated by scatter- 
ing techniques. The idea at the heart of TRC is that the 
temporal evolution of the speckle pattern generated by 
the scattered light will be very different for homogeneous 
vs. heterogeneous dynamics. For homogeneous dynam- 
ics, we expect the speckle images to change smoothly 
in time. By contrast, for heterogeneous dynamics the 
speckle pattern is expected to evolve discontinuously, be- 
cause the rate of change of the sample configuration will 
not be constant but rather fluctuate with time. Exper- 
imentally, the speckle images are recorded by a CCD 
camera and their evolution is quantifled by introducing 
c/(i, r), the degree of correlation between images taken 
at time t and t + t {a rigorous definition will be given 
in sec. Inspection of the i-dependence of c/ at fixed 
r allows temporally heterogeneous dynamics to be dis- 
criminated from homogeneous dynamics. Indeed, in the 
former case a large drop (increase) of c/ is observed when- 
ever the dynamics is faster (slower) than average, while 
in the latter the degree of correlation is constant. This 
method is quite general, since it can be applied to any 
experimental configuration where a multi-element detec- 
tor can be used to record the speckle pattern generated 
by a sample illuminated by coherent radiation. Examples 
are CCD-based light scattering experiments in the single 
scattering regime, both at wide [T^ and small angle '20|, 
DWS in the transmission or backscattering 23] ge- 
ometry, and X-photon Correlation Spectroscopy (XPCS) 
at small angles In principle, the technique could 

also be extended to non-electromagnetic radiation, e.g. 
to acoustic DWS [25l |. 

Experiments on diluted suspensions of colloidal Brow- 
nian particles have shown that the degree of correlation 
exhibits some fluctuations even in the absence of dynamic 
heterogeneity As it will be shown in detail, these 

fluctuations are due to statistical noise stemming from 
the finite number of speckles in the CCD images. In 
order to exploit quantitatively TRC data it is thus nec- 
essary to separate the contribution to the fluctuations of 
c/ due to the noise from that due to dynamic heterogene- 
ity. Although in most cases it is not possible to directly 
correct the TRC time series for the noise, we will show 
that it is possible to correct statistical quantities derived 
from the TRC data and used to quantify the fluctuations 
of c/. We focus in particular on three statistical objects: 
the time variance of ci{t, r), the probability distribution 
function (PDF) of cj{t,T) for a flxed r, and the time 
autocorrelation of the ci{t,T) trace itself. 

The variance of c/, cr^^(r), is the lowest moment of 
the data that provides information on the fluctuations. 
It corresponds to the so-called dynamical susceptibility 
X4 studied in many simulation and theoretical works on 
glassy systems 0[ [S^, H^. In a typical simulation, 
Xi is the variance of the intermediate scattering function, 
or a similar correlation function describing the system's 



change in conflguration, which is obtained from several 
independent runs. Similarly, a^^ quantifles the fluctua- 
tions of the intensity correlation function as the system 
evolves through statistically independent conflgurations. 
Importantly, a^^ allows one to relate temporal dynamic 
heterogeneity to spatial correlations of the dynamics. In 
fact, it can be shown that X4 is the volume integral of 
the spatial correlation of the local dynamics (s^l . There- 
fore, large values of a^^ will be indicative of long-range 
correlations of the dynamics. Intuitively, one can expect 
the variance of the fluctuations of the dynamics to scale 
as the inverse number of "dynamically independent" re- 
gions in the scattering volume, and thus to increase as 
the spatial range of the correlation of the dynamics in- 
creases. Recent TRC experiments on a shaving cream 
foam support this simple picture T^] . 

The PDF of the TRC signal is the most complete sta- 
tistical characterization of the dispersion of the data. 
Any deviation from a Gaussian shape immediately hints 
to heterogeneous dynamics, as suggested by experi- 
ments on a variety of systems, including colloidal gels 
[21t. |31j ^] , concentrated colloidal suspensions jS^I and 
surfactant phases [s^], foams [s^, and granular materi- 
als [S^ . Remarkably, the shape of the PDF of c/ is often 
strongly reminiscent of that obtained for similar quanti- 
ties in theoretical and simulation work on other glassy 
systems. An example is provided by simulations of spin 
glasses, where the fluctuations of the correlation func- 
tion of the spin orientation are distributed according to a 
generalized Gumbel PDF [s^l , a probability distribution 
characterized by an exponential tail strikingly similar to 
those reported for c/ in refs. "27, 31, 35,|3l|3l|33 or 
shown in this paper (see flgs. ISl and ITT)) . Similar distri- 
butions are also obtained in a variety of numerical and 
analytical investigations of systems with heterogeneous 
dynamics, both at equilibrium and out-of-equilibrium 
HE HE ^ ■ Indeed, it has been proposed that the 
Gumbel distribution arises as a universal PDF for vari- 
ous quantities measured in systems with extended spatial 
and/or temporal correlations [l^l. Clearly, in order to 
compare in detail and quantitatively the PDF measured 
in TRC experiments to those obtained analytically or by 
simulations it is necessary to correct the experimental 
data for the contribution of the measurement noise. 

The variance and the PDF of c/ describe the disper- 
sion of the data, but are insensitive to the way the fluc- 
tuations are distributed in time. By contrast, the time 
autocorrelation of the TRC signal, which was introduced 
in ref. [33| and which we shall term "second correla- 
tion" , provides information on the temporal organization 
of the fluctuations of the dynamics, and sheds light on 
the rate and life time of rearrangement events. The sec- 
ond correlation introduced here is similar to the 4-th or- 
der intensity correlation function proposed by Lemieux 
and Durian l43l| and to the multi-time correlation 
functions measured in nuclear magnetic resonance exper- 
iments probing dynamical heterogeneity near the glass 
transition ^4^. Moreover, we note that the second cor- 
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relation is the analogous, in the time domain, of the so- 
called second spectrum originally introduced to investi- 
gate no n-g aussian fluctuations in the dynamics of spin 
glasses [43 ■ 

In this paper, we first describe how to optimize a TRC 
measurement by correcting the CCD data for the dark 
noise, due to the electronic noise of the CCD, and for the 
uneven illumination of the detector (sec. We then 
turn to the temporal fluctuations of c/ and analyze sep- 
arately the contribution of the measurement noise, due 
to the finite number of pixels, (sec. IIIip and that due to 
heterogeneous dynamics fsec. lIV|) . We illustrate our anal- 
ysis by showing TRC data obtained mainly from a dilute 
suspension of Brownian particles and a shaving cream 
foam, two model systems that exhibit homogeneous and 
heterogeneous dynamics, respectively. We then intro- 
duce a method for correcting the experimentally mea- 
sured variance and PDF of c/ for the noise contribution. 
The method consist of an extrapolation scheme based on 
the observation that the measurement noise vanishes in 
the limit N ^ 00 (TV being the number of CCD pixels 
used to calculate c/), whereas the fluctuations due to dy- 
namical heterogeneities are independent of N. In sec. 
we introduce the second correlation and provide a work- 
ing formula for correcting it for the noise contribution. 
Possible artifacts that may affect the fluctuations of c/ 
and lead to spurious contributions to its PDF and to the 
second correlation are discussed in sec. I VII together with 
methods to detect or correct them. In the concluding sec- 
tion, we briefly discuss the connections between TRC and 
other techniques that measure dynamical heterogeneities 
in soft glassy systems, focussing on the advantages and 
the limitations of the various methods. 



II. TIME RESOLVED CORRELATION: 
MEASURING THE TIME-DEPENDENT 
DEGREE OF CORRELATION ci{t,T) 

In a TRC experiment a CCD camera is used to 
record, at constant rate, the speckle pattern of the light 
scattered by the sample [l^ . The CCD images are stored 
on the hard disk of a personal computer (PC) for later 
processing. The degree of correlation, c/, between speckle 
patterns at times t and t + r is then calculated according 
to 



c/(t,r) = 



G2(t,r) 

I{t)I{t + 7 



- 1 



(1) 



where G2(i,T) = {Ip{t)Ip{t + t))p and I{t) = {Ip{t))p, 
with Ip{t) the intensity measured at time t by the p-th 
CCD pixel and (• • ■)p an average over N pixels. Note 
that the normalization factor in Eq. (Q, I{t)I{t + r), 
allows to cancel out exactly any fluctuations due to 
changes in the laser power. For stationary dynamics, 
the intensity correlation function usually measured in a 
light scattering experiment, g2{T), may be obtained from 
52(-r) - 1 = ci{t,T). 



In Eq. the intensity value at each pixel is required. 
In practice, however, the images recorded by the CCD are 
affected by a pixel- and time-dependent dark (or elec- 
tronic) noise, Dp{t), and, possibly, by the non- uniform 
illumination of the detector, depending on the setup op- 
tics. To account for these effects, we write the raw signal 
for pixel p, Sp, as Sp{t) = Ip{t)bp/ {bp)p + Dp{t), where 
bp is a time- independent, spatially slowly- varying func- 
tion that accounts for non-uniform illumination and {bp)p 
is introduced so that the factor multiplying Ip(t) has 
unit average. Prior to each measurement, we collect 100 
dark images by covering the CCD detector with a black 
cap. These dark images are used to calculate the time- 
averaged dark noise Dp. The raw signal is corrected by 
subtracting, pixel-by-pixel, the average dark noise. To es- 
timate bp, we average over time the dark-noise-corrected 
CCD signal : bp = Sp{t) — Dp. For experiments whose 
duration, T^xp, is much longer than the relaxation time of 
92{t), To, this procedure averages out the spatial fluctua- 
tions associated with the speckles and leads to a smooth 
function bp, since for each pixel the intensity fully fluctu- 
ates many times and its time average is a good estimator 
of the mean intensity at a given location of the CCD 
detector. By contrast, when Tcxp ^ tq, bp keeps some 
memory of the speckled appearance of the instantaneous 
intensity distribution. In this case, we further smooth bp 
by averaging it over a few adjacent pixels. The desired 
intensity to be used in Eq. ^ is obtained from 



Ip^i) = [Sp{t)- D,;\{bp)p/bp. 



(2) 



Note that Ip{t) is still affected by an instanteneous noise 
tp{t) = Dp{t) — Dp, since only the average dark noise 
could be subtracted. By definition, = 0, while the 
standard deviation of e is typically of the order of 
1/100 of the saturation level for an 8-bit CCD camera. 
When all the relevant time scales of the dynamics are 
much longer than the inverse CCD rate, ep(t) may be 
considerably reduced by averaging Ip{t) over a few CCD 
images before applying Eq. In the following, we will 
disregard ep{t), unless when explicitly stated. 

We show in fig. how the different corrections dis- 
cussed above affect the intensity autocorrelation function 
measured for a dilute suspension of monodisperse Brow- 
nian particles. Data are collected in the single scattering 
geometry, at a scattering angle 9 — 45 deg. The particles 
are polystyrene spheres of radius a — 265 nm, suspended 
at a volume fraction (j) — 3.7 x 10~^ in almost pure glyc- 
erine cooled at 15 °C in order to make the dynamics slow 
enough to match the limited acquisition rate of the CCD, 
which is 2 Hz in this experiment. Both the dark noise 
and the non-uniform illumination result in a spurious in- 
crease of the base line of the autocorrelation function 
(squares and circles in fig. and lead to a change of 
its r intercept (see inset). However, the relaxation 
time obtained from the small r behavior of the intensity 
autocorrelation function is essentially unaffected by the 
dark noise and the non-uniform illumination. When the 
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FIG. 1: Semi-logarithmic plot of the intensity time autocor- 
relation function measured for a suspension of monodisperse 
Brownian particles in the single scattering regime {0 — 45 
deg). Data are averaged over A'' = 151800 pixels and over 
Tcxp ~ 10800 sec. Squares: intensity autocorrelation function 
calculated from the CCD raw signal, Sp{t); circles; the same 
data after correcting for the dark noise only; triangles: after 
correcting for both the dark noise and the uneven detector 
illumination. The line is a fit to the corrected data by a sin- 
gle exponential decay, the functional form of 32 — 1 predicted 
for monodisperse Brownian particles, yielding a characteristic 
time To = 5.8 sec. Inset: zoom of the small r behavior of the 
intensity autocorrelation function. 



CCD signal is corrected according to Eq. ||2Jl (up trian- 
gles) , a single exponential decay is observed, as predicted 
for monodisperse Brownian particles. For the corrected 
data, the base line is limited only by the dark noise €p{t): 
its value is of the order of 2 x 10~^ and is comparable to 
that obtained in traditional light scattering setups, which 
use a photomultiplier tube or an avalanche photodiode as 
a detector. 



III. TEMPORAL FLUCTUATIONS OF ci{t,T): 
THE NOISE CONTRIBUTION 



When disregarding the fluctuations of ep{t), the tem- 
poral fluctuations of the degree of correlation c/(t, r) at a 
fixed lag r have only two independent sources: the statis- 
tical noise due to the finite number of speckles probed in 
the experiment and the intrinsic fluctuations of the sam- 
ple dynamics. The first contribution is always present: 
we shall refer to it as to the "measurement noise" , not to 
be confused with the dark noise discussed in Sec.|n] The 
second contribution, on the contrary, is present only if 
the dynamics is temporally heterogeneous and thus rep- 
resents the physically valuable information that we aim 
to extract from the fiuctuations of c/. To highlight the 
two different contributions, we rewrite Eq. as 

ci{t,T) = g2(ai{t), ...,afe(t);r) - 1 + n(t, r), (3) 



where n{t, r) is the measurement noise, with n{t, r) = 0, 
and g2(ai(i), afe(i); t) — 1 is the pixel-averaged two- 
time intensity correlation function that would be mea- 
sured in the absence of any noise, i.e. if c/ was averaged 
over an infinite number of speckles. ai{t), ak{t) are 
parameters that fluctuate with time if the dynamics are 
heterogeneous, but are constant for homogeneous dynam- 
ics. 

To illustrate this point, let us consider as an example 
of homogeneous dynamics the single scattering measure- 
ment of the dynamics of monodisperse Brownian parti- 
cles at a scattering vector q. If the temperature is care- 
fully controlled, the diffusion coefficient D of the particles 
does not evolve with time and 52(ai(i), afc(i); r) — 1 = 
fli exp(— a2T), with ai = (3 and 02 — 2q^D constant |l,l| . 
In this case, the fluctuations of c/ are due only to the 
noise n{t,T). By contrast, the spatially and temporally 
localized bubble rearrangements in a shaving cream foam 
provide a simple example of heterogeneous dynamics. For 
a foam, the average intensity correlation function mea- 
sured in a DWS experiment in the transmission geom- 
etry has a shape very close to that for Brownian par- 
ticles in single scattering: 'g2{T) — 1 = exp[— (27t)^]. 
with 71 « 0.9 and 7 = Rr^L^/l*"^ ^ Here, R is 

the average bubble rearrangement rate per unit time and 
unit volume, is the typical volume that is rearranged 
by a single event, L is the sample thickness, and I* is 
the photon transport mean free path. In contrast with 
the Brownian suspension, however, the degree of correla- 
tion measured for a foam fluctuates not only because of 
the noise, but also because the instantaneous rearrange- 
ment rate continuously changes due to the intermittent 
nature of the bubble dynamics 0,113. Hence, for the 
foam g2iai(t),...,ak{t);T) - 1 = oi exp[-(a2(t)T)''3(*)], 
with ai — (3 constant while a2(i) = R{t)r^ L'^ /l*'^ and 
0.3 (^) = f^{t) fluctuate with time. 

In view of the correction scheme described in sees. 
IIVI and it is useful to first analyze the contribution 
to the fluctuations of c/ due to the noise. We assume 
that the sample dynamics be temporally homogeneous 
and stationary, as for the dilute suspension of Brown- 
ian particles discussed above. In this case, the param- 
eters ai,...,afe in Eq. Q arc constant and the fluctua- 
tions of cj are due only to the noise term, n{t,T). Since 
ci{t,T) is averaged over a large number of pixels (typ- 
ically N > 10000), its temporal fluctuations at a fixed 
T are expected to be Gaussian distributed, because of 
the central limit theorem (see Appendix B) . Accordingly, 
only c/ and the variance (j'^{t) of the noise n{t,T) are 
needed to obtain the full probability distribution of c/: 
Uci) = (2^a2)-i/2exp[ 

geneous dynamics, cr^(T) = 
To calculate , we recall that the variance of a quantity 
y that depends on the m variables xi, Xm is given by 



Pr,{ci) = (2^(T2)-i/2exp[-(c/ - c7)2/2ct2]. For homo- 

XT)=Cl{t,TY-Cl{t,Tf 



m r 

E 

i=l 



dy 
dxi 



E 



dy dy 
dxi dxi 



5 



where at. = xj 



is the variance of Xi and ax x 



the covariance between Xi and Xj {i 7^ j). 
The first sum accounts for the sensitivity of y to the 
fluctuations of the independent variables Xi, Xm, while 
the second sum accounts for any correlations between the 
Xi's. If distinct x^'s are uncorrelated, for 
i j and the second sum vanishes. 

By applying Eq. |@J to the definition of c/, Eq. we 
find 



<(r) = l// <(r) + 2G2(r) // af 



r2 ,-r6 



'-r5 



+ 2G2(t) // fT/,j(T)-4G2(r)// aa.Ar 



(5) 



where we have introduced the notation J{t) = I{t+T). In 
writing Eq. we have used I = J and 0-02,/ — <^G2,J^ 
because the scattered light was assumed to be stationary. 

The physical origin of the fluctuations of / and G2, 
quantified by aj and cr^^, as well as that of the corre- 
lation between / and J and between / and G2, quan- 
tified by (T/^j and aG2,ii is the finite number of pixels 
over which the instantaneous intensity and the inten- 
sity correlation are averaged. To illustrate this point, let 
us consider as an example I{t). As the sample evolves 
through different configurations, the speckle pattern fluc- 
tuates with a characteristic time tq. Since the instanta- 
neous pixel-averaged intensity is calculated for a finite 
set of speckles, different speckle patterns yield slightly 
different values of /(t). The larger the number of the 
sampled speckles, the closer will I{t) be to the "true" 
value of the average scattered intensity, whose estimator 
is /, and thus the smaller will be (t|. Indeed, we show 
in Appendix A that erf ~ N~'^, as expected from the 
central limit theorem . Moreover, one expects the in- 
stantaneous pixel-averaged intensity at time t to be cor- 
related to the same quantity measured at time t + at 
least for r < tq, because it takes a few tq for the speckle 
pattern to be completely renewed. Therefore, the covari- 
ance term (T/,j(t) will vanish only for r tq. More 
precisely, we expect cr/,j(r) to be proportional to G2(t). 
In fact, (7/^,/(t) is precisely the un-normalized correla- 
tion function measured by using the whole CCD chip as 
a single detector, similarly to the case of a traditional 
light scattering experiment where the detector collects a 
large number of speckles. 

Similar arguments may be invoked for CTq^ (r) and 
CTGg ,/(■'') I suggesting that the variance and covariance fac- 
tors in Eq. ||SJ) scale with and depend linearly on the 
average correlation function: 



(^x„x,{t) = N [Ax,,x, + Bxi,x,Ci{T)], 



(6) 



where Xi and Xj stand for any of /, J, and G2, while 
Axi,xj and Bxi^xj are constants whose values are given in 
Appendix A. To test the linear dependence of the vari- 
ance and covariance factors on the time-averaged corre- 
lation function, we plot parametrically crQ^(r), aj^jir), 
and aG2,i{T) as a function of c7(t), for data obtained in 
the single scattering experiment on Brownian particles. 
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FIG. 2: Parametric plot of the normalized variance a nd co- 
variance terms in the r.h.s. of Eq. as a function of c/(r), 
for the same experiment as in fig. El The lines are linear fits 
to the data, suggesting the ci dependence of Eq. @. 



as shown in fig. El In all cases the data are very well fitted 

by straight lines, thus confirming the validity of Eq. JSJ. 

The T dependence of the variance of the measurement 

noise can be obtained by substituting Eq. |(HJ into Eq. (jSJ. 
2 

Using G2(t) = / (c/(t) -I- 1) (see Appendix B), one finds 

a third order polynomial dependence of the variance of 
n(i, r) on c/(t): 



1=0 



aiciir) 



(7) 



where the coefficients ai can be obtained from /, Ax^.xj, 
and Bxi^x ■ Note that this third-order polynomial de- 
pendence is due to the choice of the normalization of c/ 
(see Eq. (^). If the denominator was chosen to be / , 
as in traditional light scattering experiments, only the 
first term in the r.h.s. of Eq. |SJ| would be non-zero and 
cr^^ (r) cx CTq^ , i.e. the fluctuations would increase linearly 

with c/(t). Although the normalization we have chosen 
leads to a more complicated expression for af^ (r) , we re- 
mind that it suppresses spurious variations of c/ due to 
fluctuations of the incoming beam power and thus should 
be used in TRC experiments. The inset of fig. |21 shows 
a semi-logarithmic plot of aci{T) vs. r for the Brownian 
particles, confirming that the noise of c/ decreases with 
T, as indicated by Eq. Q. In the main plot, the same 
data are plotted parametrically as a function of c7(t) . A 
very good agreement is found between the experimental 
data and the polynomial form suggested by Eq. (0, as 
shown by the line. 

The dependence of a'^{T) is the key feature that 
will be exploited in the correction for the measurement 
noise. To test this scaling, we analyze the time series 
of speckle images recorded for the Brownian particle sus- 
pension, for which a"^^ = a^, by processing different num- 
ber of pixels. First, all pixels of each image are processed 
and cj{t,T) and its variance o'^^(t) are calculated. Each 
image is then divided into two regions of interest (ROI) 
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FIG. 3: Main plot: param etric representation of the variance 
of cj as a function of ci(t) for the same experiment as in figs.0 
and 121 (single scattering from a dilute suspension of Brownian 
particles). The line is a fit of Eq. (|7J to the data. Inset: same 
data plotted vs. t. 




FIG. 4: For the same experiment on Brownian particles as 
in the preceding figures: variance of c/ as a function of the 
inverse number of pixels over which c/ is averaged. The lines 
are linear fits to the data for A^~^ < 1.2 x 10"*, demonstrating 
that for large A'^ a'i^ oc . The inset zooms in the region 
near to the origin. 



of equal size. For each ROI, c/ and its variance are calcu- 
lated and the values of cr^^ (t) obtained for the two ROIs 
are averaged, yielding the variance of c/ when only N/2 
pixels are processed. This scheme is iterated until the 
size of each ROI is reduced to 225 pixels. Figure0]shows 
'^ci {''') a function of the inverse number of processed 
pixels, N^^, for three time delays corresponding to 0.09, 
0.87, and 8.7 times the relaxation time of ^2 ^ 1, respec- 
tively. In all cases, the data for N^^ < 1.2 x 10~* are 
very well fitted by a straight line that goes through the 
origin, as shown in the inset |6^ . This confirms that for 
temporally homogeneous dynamics a^^ oc N^'^, as indi- 
cated by Eq. {T)). Note that a deviation from this linear 
trend is observed at the largest N~-^, due to edge effects. 
In fact, the contribution of each pixel to c/ is not com- 
pletely independent from that of nearby pixels, because 
the intensity of the speckle pattern is spatially correlated 
over a distance of a few pixels. Pixels far from the edges 
of a ROI have more nearby pixels than those on the edges; 
accordingly, the statistically independent contribution to 
c/ carried by a bulk pixel is less than that of an edge 
pixel. When reducing the size of the processed ROI, the 
weight of edge pixels relative to bulk pixels increases and 
corrections to the scaling become increasingly ap- 
parent. These corrections are negligible for N > 8000, 
as seen in fig. 0] We find that a similar scaling is 
obtained at all time delays (data not shown). Indeed, the 
slopes of the straight line fits to tr^^ vs. provide an 

estimate of the proportionality coefficient '^i^QCtici{T) 
that agrees within 1% with the value directly obtained 
when calculating a^^ by processing the maximum number 
of available pixels. 

For temporally homogeneous dynamics, the statistics 
of the variables 6*2 (i,T) and I{t) in Eq. Q is Gaussian, 




0.48 0.52 0.20 0.24 0.00 



FIG. 5: Symbols: PDF of ci{t,T), Pcj{ci), for the same sus- 
pension of Brownian particles as for the preceding figures. 
From left to right, the normalized delay t/tq is 0.09, 0.87, 
and 8.7. In the three panels the cj axis spans an interval 
of equal width (0.04). The lines are Gaussian distributions 
whose mean and standard deviations are directly obtained 
from the ci time series, without any fitting parameter. 



because of the central limit theorem. As a consequence, 
the probability density function (PDF) of c/ is also Gaus- 
sian, as demonstrated in Appendix B. In fig. |Slwe show 
the PDF of c/ for various r for a Brownian suspension 
of particles. The symbols are the experime ntal d ata, 
while the lines are Gaussian PDFs with mean c/(t) and 
standard deviation CTc/ (t) obtained directly from the c/ 
time series, without any fitting parameters. An excellent 
agreement between the data and the theoretical shape of 
the distributions is observed at all r. 
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IV. TEMPORALLY HETEROGENEOUS 
DYNAMICS: CORRECTIONS FOR THE NOISE 
CONTRIBUTION 



In temporally heterogeneous dynamics, the fluctua- 
tions of c/ are due both to the noise and to dynamical 
heterogeneity. In this section wc propose a method for 
correcting the variance and the PDF of c/ for the noise 
contribution, so as to obtain the statistics of the fluctua- 
tions due to dynamical heterogeneity. Moreover, we show 
that in some instances the full temporal evolution of c/ 
may be corrected for the noise, thus allowing 32(^1 ''") ^ 1 
to be determined, not only its variance and PDF. 



A. Correction of the variance of ci 

In the case of dynamically heterogeneous processes, the 
parameters oi, aj, in the two-time correlation function 
52 (oi ) , . . . , flfc ) ; r) — 1 , Eq. Q , fluctuate with t. There- 
fore, an extra term, <Tg^ , contributes to the variance of c/, 
in addition to the noise term analyzed in the preceding 
section: 



Eq. lO yields, for the foam. 



92 



(8) 



l=Q 



In writing the expression above, we have assumed that 
no correlation exists between the noise n{t, r) due to the 
finite number of pixels and the fluctuations of ai, Ofc. 
Using Eq. cr^^ may be expressed as 



E 

i=l 

E 



992 

dai 

dg2 dg2 



(9) 



An example where assumes a particularly simple form 
is given by the dynamics of a shaving cream foam, result- 
ing from intermittent bubble rearrangements, as mea- 
sured in a DWS experiment. We find that the fluctua- 
tions in the instantaneous decay rate ^{t) of the correla- 
tion function are slow compared to 7, so that at any given 
time (72 — 1 is well approximated by a stretched exponen- 
tial, /3exp[— (7(t)r)^(*']. Small variations of account 
for slight changes of the decay rate on a time scale com- 
parable to 7~^. For example, if the dynamics tend to 
slow down during the measurement of 32 — 1, the initial 
decay of the correlation function will be faster than its 
final decay. Thus, the shape of §2 will be more stretched 
than the average one {fj, < Jl). Conversely, > /I if the 
dynamics accelerate during the measurement of 172 — 1- 
By taking into account the fluctuations of both 7 and fi, 



<M = (7r)2/32exp[-2(7T)^] 



[ln(7r)cr^] 




(10) 



with (7^ = (7 - 7) 



7)2 and (T^ — (/X — 71)2. Equation ifTUI) 
extends a similar expression given in ref . [T^ l , where only 
the fluctuations of 7 where taken into account. However, 
here we still neglect possible correlations between 7 and 
/i that could be described by including the second term 
of the r.h.s. of Eq. ®. 

A direct test of Eq. (|10|) is not possible, since we ex- 
perimentally only access a^^ , not a^^ . Therefore, for the 
foam as well as for the general case of temporally hetero- 
geneous dynamics it is desirable to subtract the trivial 
contribution of the measurement noise from the experi- 
mentally measured cr^^ , in order to obtain the physically 
relevant variance . We have tested two different ap- 
proaches. In the first method, the linear dependence of 
(Jxi,xj on cj shown in fig. [21 is used to derive formulas 
for these quantities that are independent of the instan- 
taneous dynamics and thus of the homogeneous vs. het- 
erogeneous nature of the dynamics. These formulas are 
presented in Appendix A. Unfortunately, although they 
provide separately fairly good estimates of the various 
terms ,x , when combined using eq. |31 to evaluate ti^ 
uncertainties add up leading to errors of the order of 
30 — 40%, as tested on the data for the single scattering 
experiment on the Brownian particles. 

In the following, we describe in detail a second method 
that has proven to be highly effective. The key point is 
to recognize that, contrary to the noise term, the fluctu- 
ations of 52(01, Cfc; r) do not depend on the number of 
pixels over which c/ is averaged. This is because in the 
far field geometry of the scattering experiments described 
in this work, each CCD pixel collects light scattered by 
the whole illuminated sample. Thus, any spatial or tem- 
poral heterogeneity of the dynamics affects in the same 
way the signal measured by each pixel. The different 
pixel-number dependence of the noise and the fluctua- 
tions (flrst and second term in the r.h.s. of Eq. |(HJ), re- 
spectively) suggests a way to discriminate between these 
two contributions. We analyze the speckle images by 
processing different number of pixels, as described for the 
Brownian suspension in sec. lIIII and plot a^^ (r) as a func- 
tion of A^^^, as shown in fig.lHl As indicated by Eq. ||SJ|, 

the slope of a linear fit to the data yields X];=o ^^'^/(t) , 
while the intercept at A^~^ = is a^^ (r), the desired vari- 
ance of the correlation function due to dynamical hetero- 
geneity. Thus, the representation of figs. 01 and [HI allows 
one to extrapolate tr^^ to A'^ = 00, where the measure- 
ment noise vanishes. As seen in the inset of fig. O for 
T7 ^ 1 or T7 ^ 1 the intercept of the linear flt is very 
close to zero, indicating that at these delays the fluctu- 
ations of cj are mainly determined by the measurement 
noise. By contrast, at intermediate delays the intercept 
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clearly departs from zero, revealing the intermittent na- 
ture of the dynamics of the foam. 

In fig. 13 we plot the r dependence of cr^^ measured 
for the foam, together with the noise contribution, a'^ = 

N~^'^i^QaiCi{T) , and that of the intrinsic fluctua- 
tions, (Tg^- The noise contribution is extracted from the 
slope of cr^^ vs. N~^, while the variance of the intrinsic 
fluctuations is given by the N^-^ limit of cr^^ vs. N~^. 
At all T we find cr^^ = + a'^^ within 0.8%, thus con- 
firming that our analysis allows to correctly account for 
the two contributions to the fluctuations of c/. Note that, 
while the shape of the time-averaged correlation function 
^2 — 1 is almost the same for the foam and the Brownian 
particles (a slightly stretched exponential for the former 
and a simple exponential decay for the latter), the r de- 
pendence of the fluctuations of c/ is very different (com- 
pare the inset of fig.Oto fig.©, thus allowing temporally 
heterogeneous dynamics to be unambiguously detected. 
For the foam, correcting a^^ for the noise contribution 
is especially important at time delays far from the mean 
relaxation time, because the intrinsic fluctuations die off 
for T ^ 0, when virtually no rearrangement had a chance 
to occur, and for r ^ oo, when so many rearrangement 
events occurred that the statistical fluctuations of their 
number are negligible. By contrast, we recall that cr^ 
remains finite at all r. Once corrected for the noise con- 
tribution, cr^^(T) is very well described by Eq. (^UJ), as 
shown by the line in fig. [7| Interestingly, the fluctuations 
are maximal on the time scale of the mean relaxation 
time, a general feature found also in the dynamic suscep- 
tibility Xi measured in simulations. Intuitively, this can 
be explained by recognizing that the correlation function 
is most sensitive to a change of the instantaneous relax- 
ation time for t « tq. Finally, note that when comparing 
the absolute values of (7^^ (r) and Xi one should keep in 
mind that the latter is usually defined as the variance 
of the correlation function multiplied by the number M 
of particles in the system. Therefore, for homogeneous 
dynamics X4 ~ 1, since the variance of the number fluc- 
tuations is of order 1/M, while X4 > 1 for heterogeneous 
dynamics. In the case of the foam shown in fig. tak- 
ing M as the number of bubbles in the scattering volume 
leads to Ma'^^{T) ^ 1000, indicating strongly heteroge- 
neous dynamics. 




0.0 1.0x10"^ 2.0x10'^ 

n' 

FIG. 6: Variance of ci{t,T) for a foam, as a function of the 
inverse number of pixels over which ci is averaged, for three 
time delays r. The lines are linear fits to the data for A'^"^ < 
5.5x 10~*. For t'j 1, the intercept of the fit strongly departs 
from zero, indicating temporarily heterogeneous dynamics as 
discussed in the text. Note that the fit yields a non-zero 
intercept also at small delay, as shown in the inset that zooms 
into the small a^^, large A'^ region. Data were collected in the 
DWS transmission geometry for a duration Texp ~ 160 sec, 
242 times longer than the relaxation time of ci, = 0.33 
sec. 




T(sec) 

FIG. 7: Variance of ci(t,T) for a foam, as a function of r. 
Open squares: variance of the raw data (c/); open circles: 
noise contribution; solid triangles: a^.^{T), the contribution 
due to the fluctuations of the number of rearrangement events 
per unit time. The line is a fit to (jg^ ij) using Eq. Data 
were collected in the DWS transmission geometry. 



B. Correction of the PDF of ci 

In the absence of intrinsic dynamical fluctuations, 
knowledge of the average degree of correlation, cj, and 
of the noise variance, cr^, is sufficient to fully determine 
the PDF of c/ at fixed lag, because c/ is a Gaussian 
variable. Heterogeneous dynamical processes, on the 
contrary, lead in general to non-Gaussian distributions 
of Ci{t,T) 1^ 31, 34], whose shape depends on both 
the dynamical process and the time delay r. Because 
ci{t, t) is the sum of two uncorrelated random variables. 



.g2(ai(<), afe(t); r) — 1 and the noise n{t,T), the PDF 
of c/ is the convolution of the probability distribution of 
52 (oi (i) , ■ • ■ , Ofe (t) ; r) — 1 with that of n(t, r) ■ 

Pe,(c/) = Pg,-i(52-l)®Pn(n), (11) 

where Px{x) denotes the PDF of x and (/ (8) g){x) — 
J dx' f{x')g{x — x') is the convolution product. In order 
to recover the physically relevant PDF of 52 ~ 1 from the 
measured Pc/(c/), one may use standard Fourier trans- 
form techniques to deconvolute the experimental data. 
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using Pn{n) as the response function Pg^-i = 

T^\T{Pci)l T{Pnf\, where T and indicate direct 
and inverse Fourier transform, respectively. Unfortu- 
nately, this procedure is very sensitive to noise in the 
data and leads typically to unstable solutions exhibiting 
wide oscillations. Instead, we use a technique similar to 
the indirect Fourier transformation (IFT) method used 
to process static scatte ring data. Details on the IFT 
method can be found in |49l l50|. here we simply describe 
the main steps of our implementation. We first assume 
that the unknown PDF oi gi — \ at fixed lag r may be 
written as the linear superposition of a set of suitable 
functions (/)/c: 

M 

=^Mfc'/'fc(52-l). (12) 

k=\ 

It is convenient to choose (f>kix) — rect(a;— x^; 2A), where 

rect(.-x.;2A).{j ; ^ (13) 

is a square pulse of width 2 A centered on Xk- 2 A is 
taken to be the width of the bins used to calculate the 
PDF of c/, i.e. the separation between the c/ coordinates 
of the M experimental Pc, data. Because the convolution 
is a linear transformation, Eqs. Hll|l and l|12|) yield the 
following guess for the PDF of c/: 

M 

Pgncss{ci) = ^ ^ferect(c/ - c/,fc; A) (g) P„(c/) 
fc=i 

M I — 

i=k ' 

with 

x^^^al (ci - ci^k ± A) . (15) 

Here erf(a;) = 27r~^/^ exp(— u^)(iu is the error func- 
tion and ci^k is the center of the fc-th bin 
used to calculate the experimental PDF. In writing 
the last hue of Eq. ((CTl we have used P„(c/) — 
(27rcr^)~^/^ exp[— Cj/2ct,^] and calculated exphcitly the 
convolution product. We note that the width of the rect 
functions is typically much smaller than that of the Gaus- 
sian noise P„, so that the difference of the erf functions 
in square brackets is very close to a Gaussian. We stress 
that the only unknowns in Eqs. H14|l and IjlSI) are the 
coefficients ^k, since the variance of the noise cr^ can 
be obtained directly from the experimental c/ using the 
extrapolation scheme described in the preceding subsec- 
tion. 

In principle, the coefficients fik can be determined by 
fitting Eq. H14|l to the experimentally measured PDF of 
c/. Once the /i^'s are known, the desired PDF of g2 — 
1 can be obtained by using Eq. (|12|l . In practice, two 



issues must be addressed when determining the set of 
/ifc. Firstly, noise in the experimental Pc/ can make the 
fitting procedure unstable. It is therefore convenient to 
smooth the experimental Pcj before fitting it by Eq. I|14|l . 
To this end, Pcj is approximated by a smooth curve, Psm, 
obtained by fitting the data by a suitable function. We 
find that in most cases 

Psm(c/) = y{Q'^-i-i'^i-'^o)+i^2ici-co)^-U3Cxp[}^4,(ci-co)] f^-^Qj 

fits well the whole experimental PDF, although fitting 
Pcj piecewise may be sometimes necessary. Note that the 
GumbcU PDF corresponds to the case vi = 1/4 = a > 0, 
V2 = 0, and = 1 [s^l- The coefficients fik are then 
found by fitting Pguoss(c/) to Psm, rather than directly 
to Pcj- Secondly, the PDF of 52 ~ 1 determined by in- 
serting the /ife's thus obtained in Eq. (|14|) often exhibits 
large oscillations. This is because the erf functions in 
Eq. 1141 1 do not form an orthogonal basis, as discussed in 
ref. • This problem can be solved by adding a stabi- 
lization condition when fitting Pgm by Pguess- We follow 
ref. and determine the set of Hk by minimizing the 
following expression: 

L M-l 

Y^[Psm{ci,l) - Pgucss(c/j)]^ + ^ H ^^'"'+^ ^ Mfc)^-(17) 
1 = 1 k=l 

The first term in the above expression corresponds to the 
usual sum of squared deviations between the fitting func- 
tion (Pguoss) and the data (Psm)- The second term assigns 
a cost to any large variation between successive /Jk and 
thus tends to suppress all fast oscillation of Pguess- The 
relative weight of the two terms is controlled by the La- 
grange multiplier A, whose optimum value is determined 
as described in ref. 's^. 

The top panel of fig.|Slshows the PDF of cj measured 
in the DWS TRC experiment on foam, for t = 0.02 sec 
(open circles). The solid line is the smoothed PDF, Psm, 
obtained by fitting the data to Eq. (|16l) . The dotted line 
is the Gaussian PDF of the noise, P„, whose width is (t„ 
(for display purposes, P„ has been centered on cJ, rather 
than on 0). In the bottom panel, the thick line is the PDF 
corrected for the noise contribution, Pg^-i- Note that 
the right wing of the corrected PDF drops much more 
abruptly than that of the raw data (for comparison, the 
uncorrected and the smoothed PDF are also plotted in 
the bottom panel). The left wing, on the contrary, is al- 
most unaffected by the correction. This is a consequence 
of the nearly exponential behavior of the left wing: in- 
deed, it can be shown that the convolution of an exponen- 
tial function with a Gaussian is again exponential, with 
the same growth rate. We test how close the raw data 
and the corrected PDF are to a Gumbel distribution by 
fitting the data to the expression of Eq. IjlGI) , with V2 = 
and = v^Vi, corresponding to a generalized Gumbel 
PDF [33. For the raw data, we find 1^3 — 1.04, very close 
to 1/3 = 1, the value for a Gumbel PDF. By contrast, for 
the corrected data = 0.42, showing that the right wing 
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FIG. 8: Top panel; PDF of ci{t, r) for a foam, for r = 0.2 sec 
(open circles); data were collected in the DWS transmission 
geometry. The continuous line is the smoothed distribution, 
Psm, obtained from a fit to Eq. 1161 . The dotted line is the 
Gaussian PDF of the noise. Bottom panel: the corrected 
distribution, Pg^-i, (thick line), together with the raw data 
and the smoothed PDF (same symbols as in the top panel). 
Note the much steeper decay of the right wing of the corrected 
PDF. The inset shows the data around the peak in a linear 
axis graph. 

of the corrected PDF strongly departs from both a Gum- 
bel distribution and the "universal" PDF of ref. 'Zl^, for 
which 1^3 — tt/2 « 1.57. A more detailed investigation of 
the shape of the PDF of c/ for various delays r and dif- 
ferent systems will be presented elsewhere: here we just 
stress the importance of the noise correction in view of 
any quantitative comparison. 

C. Direct correction of c/(t,r) for r << tq 

The methods developed in subsecs. llV Al and IIV Bl al- 
low one to calculate the variance and the PDF of 32 — 1, 
but not to correct directly the time-dependent degree of 
correlation ci{t,T). Here we show that such a correction 
is possible, i.e. that the noise- free g2{t,T) — 1 at fixed 
r may be retrieved as a function of t, provided that the 
following assumptions are fuUfilled: i) the dynamics is 
homogeneous on a time scale comparable to the CCD 
exposure time; ii) t << tq, where tq is the average re- 
laxation time of 52 ^ 1- 

We first observe that ci{t,T = 0) measures the so- 
called contrast of the speckle pattern, or coherence factor 
p. The latter is determined only by the speckle-to-pixel 
size ratio, which is a time-independent quantity, and by 



the blurring due to the fluctuations of the speckle during 
the time the CCD chip is exposed 0|. If i) is fulfilled, 
the amount of blurring is constant over time and hence 
c/(i, 0) fluctuates only because of the noise |6^. There- 
fore, n{t, 0) can be directly obtained from the experimen- 
tally measured c/(t,0): 

n{t,0) = ci{t,0) - ci{t,0) (18) 

If in addition also ii) is fulfilled, n{t,T) « n{t,0), be- 
cause the noise evolves on the same time scale as the 
speckle pattern, tq. Indeed, by analyzing TRC data 
for Brownian particles, we have shown in ref. f3| that 
the noise is highly correlated for r << tq. This sug- 
gests that ci{t, t) may be directly corrected according to 
52(i, t) - 1 = c/(t, r) - n{t, r) « c/(t, r) - n{t, 0), with 
n{t, 0) obtained via Eq. (|18|l . This approximation may be 
further refined by taking n{t, r) k, [n{t, 0) -I- n{t + r, 0)]/2 
rather than n{t, t) « n{t, 0) and by scaling this estimate 
of the noise so that its standard deviation matches the 
actual standard deviation of n{t, r): 

g2{t,T)-l = ci{t,T) -[n{t,0) +ji{t + T,0)]/2 

X a„(r)/a„(0). (19) 

Here, the standard deviation of the noise of ci{t,T), 
a„ {T), is obtained by applying the extrapolation scheme 
described in subsec HV Al while n{t, 0) and n{t + T, 0) are 
given by Eq. H18|l and cr„(0) is directly calculated from 
n{t,0). 

We test Eq. ((T^ on TRC data taken for the dilute sus- 
pension of Brownian particles, for which the fluctuations 
of ci are due only to the measurement noise: ci{t,T) = 
•n{t, t) -I- ci{t, t). The top panel of fig. El shows a portion 
of c/(i,0) and c/(i,Ti), with ri = 0.5 sec << tq = 5.8 
sec. Clearly, the two traces are highly correlated, as ex- 
pected if n(t, Ti) K, n{t, 0). The dotted line is g2{t, Ti) — 1 
calculated according to Eq. fl^ . Notice that g2{t, n) — 1 
is almost constant, as expected for temporally homoge- 
neous dynamics, thus demonstrating the effectiveness of 
the correction. The bottom panel shows the PDF of 
c/(t, Ti) and g2{t,Ti) — 1 calculated over the whole du- 
ration of the experiment. By correcting c/(i,ri) for the 
noise, its variance is reduced by almost a factor of 100 
(from 3.0 x 10"^ to 3.9 x 10"^). The residual fluctuations 
of g2{t, Ti) are most likely due to the noise ep{t) discussed 
in reference to Eq. ^ , whose variance we estimate to be 
of the order of 3.9 x lO"'^ [66l |. 

D. Influence of the duration of the experiment 

The analysis of the fluctuations of c/ presented in sub- 
secs. IIV Al and llV Bl was developed under the assumption 
that the dynamics be stationary, so that data could be 
collected over a period, Toxp, much longer than the aver- 
age relaxation time of the intensity correlation function, 
Tq. Dynamical heterogeneity, however, appear to be more 
prominent for systems close to jamming or quenched in 
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FIG. 9: Top panel: ci{t,0) and ci{t,Ti = 0.5 sec) measured 
for a diluted suspension of Brownian particles, for which the 
relaxation time of 32 — 1 is tq — 5.8. Data are plotted over 
a short time window to show the strong correlation between 
the two traces. The dashed blue line is (72(i, ti) — 1 obtained 
by correcting ci{t,Ti) for the noise contribution according 
to Eq. 11911 : the fluctuations due to the noise are almost 
completely removed. Bottom panel: PDF of ci(t,Ti) (open 
squares) and of g2{t,Ti) — 1 (solid circles), obtained over the 
full duration of the experiment (T^xp = 10800 sec). The PDFs 
are labelled by their variance. 



an out-of-equilibrium state. For these systems, meeting 
the condition Toxp 3> tq is often impossible, since the 
sample may be aging, leading to non-stationary dynam- 
ics, or because, even if the dynamics is stationary, the 
relaxation time may be as large as several tens of hours. 
It is therefore important to address the issue of the in- 
fluence of the experiment duration on the measured . 

Let us first consider the simpler case of homogeneous 
dynamics. We divide the time series of ci(t,T) obtained 
for the Brownian particles into non-overlapping segments 
of duration T < Texp. We denote by cr^^ ^(r) the mean 
value of the variance of ci{t,T) calculated for each seg- 
ment of duration T, and plot in fig. llOl a^^ as a function 
of T normalized by the relaxation time tq (the data refer 
to r = To, a similar behavior is observed for all r). For 
T/tq >> 1, cr^j rp is independent of T, because in this 
regime the experiment duration is long enough for the 
system (and thus the speckle pattern) to sample a suf- 
ficiently large number of different configurations. Con- 
sequently, tr^^ rp saturates to the maximum value given 
by Eq. Q, which is dictated only by r and the num- 
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FIG. 10: Variance of ci{t, t — tq) calculated over a duration T 
as a function of T normalized by the mean relaxation time tq , 
for a suspension of Brownian particles. A similar behavior is 
observed for all r. The line is the contribution to the variance 
of ci{t,T = To) of the CCD electronic noise (see the text for 
more details). 



ber of CCD pixels. Note that for T/tq close to one 
(T/to = 0.87), (T^^ rp is still significantly lower than its 
saturation value (about 35%), while it increases to 95% 
for T/tq > 20 and saturates only for T/tq > 50. As T/tq 
decreases below 1, the amplitude of the fluctuations is 
significantly reduced when decreasing T, since the sys- 
tem is not given enough time to explore significantly dif- 
ferent configurations. For T/tq ^ 1, the speckle pattern 
is essentially frozen on the time scale of the experiment 
duration: the fluctuations of c/ due to the evolution of 
the speckle pattern are thus expected to be almost com- 
pletely suppressed. For the data shown in fig. ^1 this 
regime is not quite reached, since the CCD acquisition 
rate could not be fast enough to allow for several im- 
ages to be acquired on a time scale much smaller than tq 
(the smallest lag between images is 0.5 sec = To/ 11. 6). 
In the T/tq regime, the main contribution to a^^ rp 
should come from the CCD electronic noise ep{t), whose 
variance, cr^, is shown as a line in fig. ^| The data 
shown in this figure clearly demonstrate that care must 
be taken when comparing the fluctuations measured in 
experiments whose relative duration Tqxp/to is different. 

For heterogeneous dynamics, we expect the behavior 
of (7^^ p to be qualitatively similar, although the situa- 
tion will be in general more complicated. In fact, in this 
case the relevant time scales to which Texp has to be com- 
pared are not only the mean relaxation time of 32 — 1, tq, 
but also the characteristic time of the "intrinsic" fluctu- 
ations of the dynamics, described by the variation of the 
parameters ai(t), ...afc(t) in Eq. Q. As an example, for 

the foam data shown in figs.Eland[7lTo ~ ^{t) = 0.32 
sec, while the temporal fluctuations of the decay rate oc- 
cur on a much longer time scale, Tfluct ~ 7 sec 's^l (see 
also fig. 1131 below and the associated discussion). Hence, 
in order to measure precisely the fluctuations the exper- 
iment should last more than about 20Tfluct, rather than 
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FIG. 11; a); time variation of c/ for a multilamellar vesicle gel, 
for r — 100 sec. The mean relaxation time is to = 5050 sec, 
the duration of the experiment is T = 80000 sec — 15.8 x tq. 

a) : raw data; b): the same data after applying the correction 
Eq. H19^ : c): PDF of the traces shown in a) (dotted line) and 

b) (solid line). 



just 20to. For other jammed materials, on the contrary, 
the intrinsic fluctuations may be much faster than tq, 
as demonstrated by the sudden sharp drops of c/ resuh- 
ing from intermittent rearran gem ents in a closely packed 
multilamellar vesicle system pi], shown in fig. Illb. In 
this case, at least for the smallest time lags, the analysis 
proposed in subsec. IIV Cl allows the slow fluctuations due 
to the measurement noise to be suppressed, while pre- 
serving the fast drops of c/ that contain the physically 
relevant information on the dynamical intermittency, as 
seen in fig. [TTb . Note the dramatic change of the shape 
of the PDF of c/ before and after correcting the data, as 
shown in fig. Illb . 



V. THE SECOND CORRELATION 

In the previous sections, the fluctuations of c/ were 
analyzed in terms of their variance and PDF. Additional 
information on the system dynamics can be obtained by 
studying not only the probability distribution of the fluc- 
tuations and its second moment, but also the way these 
fluctuations occur in time. We characterize the temporal 
properties of c/ at a fixed lag r by introducing the time 



Ce,(At,r) 



Cl{t,T)ci{t + At,T) 



ci 



(20) 



With this choice of the normalization, Ccj^(0,t) = 1, and 
Cci{At,T) = if ci{t,T) and ci{t + Ai, r) are uncorre- 
lated (e.g. for At — > oo). Because Ccj is the correlation 
function of a time- varying quantity -ci(t,T)- which is 
obtained itself by correlating the scattered intensity, we 
shall term it the "second correlation" , in analogy with 
the "second spectrum" first introduced in the context 
of spin glasses The second spectrum describes in 

the frequency domain the fluctuations around the mean 
value of the spectrum of a time-dependent quantity. Sim- 
ilarly, the second correlation describes — in the time 
domain — the fluctuations of the degree of correlation 
between the time-dependent system configurations. The 
second correlation is also similar to the 4-th order in- 
tensity correlation function introduced in refs. ^3.E3j 
gfj^\T) = I{0)I{T)I{t)I{t + T)/'f. Note however that 

(t) compares the scattered intensity, and thus the 
sample configuration, at four successive times, while the 
second correlation compares the change in sample con- 
figuration occurring over two time intervals of duration 
T separated by a time At. 

As for the case of the PDF discussed in sec. IIVI 
the second correlation contains contributions from 
both the physically relevant intrinsic fiuctuations of 
g2{ai{t), ...,ak{t);T) — 1 and the noise, n{t,T). We 
focus on stationary dynamics measured over a pe- 
riod Texp ^ To and assume that the fluctuations of 
g2{ai{t), ...,akit);T) — 1 and n{t,T) are uncorrelated: 
'gTn — for all At. Under these assumptions and 
by using Eq. (O, Eq. yields 
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{T)C,,{At,r)+al{T)C^{At,T) 



(21) 



where Cgj and C„ are the correlation functions of the 
fluctuations of g2 and n, respectively, defined similarly 
to Eq. (|20|l . Experiments on Brownian particles (homo- 
geneous dynamics) show that C„(At, r) ~ ci{At)/(3 for 
all r Indeed, C„ is expected to be proportional to 
c7 because the r dependence of the noise time autocor- 
relation stems from the same physical mechanism lead- 
ing to the decay of the intensity correlation function, i.e. 
the renewal of the speckle pattern over time discussed in 
sec. mil In order to extract the desired second correla- 
tion of 32 from the noise-affected 0^ , we use a method 
similar to the extrapolation technique adopted for the 
calculation of the variance of g2. At At and r fixed, Ccj 
depends on the number N of processed pixels through 
the tr^j term in Eq. (|21|l . We omit for simplicity the ex- 
plicit dependence on At and t and insert Eq. into 
Eq. 
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FIG. 12: Correction method for Cg^- Symbols: the second 
correlation Cc, as a function of the inverse number of pro- 
cessed pixels, for various At (shown by the labels, in sec) and 
for T = 1.6 sec (only the data for the largest N are shown). 
The data are obtained from a DWS measurement on a foam. 
The lines are fits to the data using Eq. I|22|l : the full range 
of the fit is shown in the figure. The outcome of the fit are 
the values of Cgj (At, r) and Cn(At,r) shown with the same 
symbols in fig. 1131 



O 
0) 



1.0 



0.5 



0.0 




0.01 



0.1 1 

(sec) 



10 



FIG. 13: Time autocorrelation function of the fluctuations of 
g2 for a foam, obtained from the second correlation Ccj (At, r) 
by applying the correcting scheme described in the text. 
Small open circles: the noise-free Cg2-i; solid line: the ex- 
perimentally measured Ccj; dashed line: the noise contribu- 
tion, Cn- The large symbols are the values of Cg^-i obtained 
from the fit to the corresponding data shown in fig. 1121 Note 
the peaks of Cg^-i at At « 7 sec, and At « 14 sec, suggest- 
ing pseudoperiodic rearrangement events, as discussed in the 
text. 



The l.h.s. of this expression can be calculated by pro- 
cessing ROIs of the speckle images of different sizes, as 
explained in sees. Illll and HVl and plotted as a function of 
. as shown in fig.ll2lfor a foam. Ther.h.s. ofEq. H22|) 
is then used as a fitting function for Cc,{N), where the 
fitting parameters are the desired and C„, while cr^^ 

and X]f=o "^'^^ obtained independently from the cor- 
rection of the variance of c/, as explained in subsec. IIV Al 
By repeating this procedure for all At and r of interest, 
the full second correlation can be corrected for the noise 
contribution. We show in fig. ^| Cg^ (At, r) for a foam 
(open circles) , together with the uncorrected data (Cc^ , 
solid line) and the noise contribution (C„, dashed line). 
Remarkably, a peak is visible in Cg^, at At = At* « 7 
sec. Note that the peak is not present in C„, whose only 
relevant time scale is that of the relaxation of cj, 7"^ — 
0.32 sec. Therefore, the peak in Cg^ must be associated 
with the intrinsic fluctuations of the dynamics due to the 
intermittent bubble rearrangements. Indeed, the peak 
indicates that the fluctuations of the instantaneous decay 
rate 7(t) are pseudoperiodic on a time scale of the order 
of At* . We are currently investigating the origin of this 
feature. 



VI. POSSIBLE ARTIFACTS 

In a TRC measurement, the focus is on the fluctu- 
ations of the degree of correlation, rather than on its 
mean value. As a consequence, TRC measurements are 
very sensitive to various sources of spurious fluctuations, 
whose effects are usually less prominent or even negligi- 



ble in traditional dynamic light scattering experiments, 
since they tend to average out. An example is provided 
in the top panel of fig. 1141 which shows c/ measured for 
the speckle pattern generated in the transmission geom- 
etry by a ground glass, a scatterer whose dynamics are 
completely frozen-in. In this case, we would expect c/ 
to fluctuate slightly only because of the electronic noise 
ep(t) discussed in sec. ^ Surprisingly, sharp drops of 
the degree of correlation are clearly visible. Because the 
drops are rare, they do not affect significantly cJ; how- 
ever, they do change significantly the PDF of c/ and the 
second correlation. Indeed, one would mistakenly take 
the dynamics to be intermittent, if the sample was not 
known. The origin of this artifact is the laser beam point- 
ing instability. Beam pointing instability is the charac- 
teristic noise of lasers that results in small fluctuations 
of the propagation direction of the output beam. Since 
the speckle pattern is centered around the propagation 
direction, any change of the incoming beam direction en- 
tails a rigid shift of the speckle image with respect to the 
CCD detector. Thus, the intensity at each pixel slightly 
changes, leading to a drop of c/. To demonstrate that 
this mechanism is indeed responsible for the sharp drops 
of c/ observed in fig. 1141 we show in the bottom panel 
the corresponding shift, Ar(t, r), between pairs of images 
taken at time t and t -I- r, measured by Particle Imaging 
Velocimetry (PIV) (s^. This method is based on spa- 
tial cross-correlation techniques and allows the rigid mo- 
tion between two images to be quantified with sub-pixel 
resolution (our adaption of PIV to speckle imaging will 
be described elsewhere). A comparison between the two 
panels of fig. El clearly shows that the anomalously large 
drops of c/ are due to larger-than-average rigid shifts of 
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FIG. 14: Top panel: ci(t,T) measured for a static scatterer, 
a ground glass, for t — 5 sec. Note the sharp drops of the 
degree of correlation that greatly exceed the noise of the data. 
Bottom panel: rigid displacement between the pair of speckle 
images used to calculate c/(t, r) above, measured by PIV. The 
sharp drops of c/ are artifacts due to larger-than-average rigid 
shifts of the speckle images. 



the speckle images. Note that shifts as small as a fraction 
of pixel, corresponding to a few microns, have a measur- 
able impact on c/. 

This example shows how sensitive to instabilities TRC 
is. Therefore, care must be taken in order to minimize 
possible instabilities and to identify any spurious fluc- 
tuations of c/. Moreover, attempts should be done to 
correct cj for artifacts. In our experience, the most com- 
mon artifact encountered in TRC experiments is simi- 
lar to that exemplified by fig. El fiuctuations or sharp 
drops of c/ due to a rigid shift of the speckle pattern, 
rather than to the characteristic "boiling" or fiickering of 
the speckle images associated with the evolution of the 
sample configuration. In addition to beam pointing in- 
stability, temperature variations are often to blame for 
spurious features in the temporal evolution of c/. 

Temperature variations may change the direction of 
propagation of light because of refractive effects, since 
the refractive index, tid, varies with temperature, T. As 
an example, let us consider a typical small-angle single 
scattering measurement on an aqueous sample contained 
in a rectangular cell, whose entrance and exit walls are 
perpendicular to the incoming beam. Light scattered at 
an angle 9s equal to, e.g., 15 deg is refracted at the water- 
to-air interface and exits at an angle given by Snell's law: 
^air = arcsin (nD sin ^s) j67j . A temperature fluctuation 



5T induces a variation (50air given by 



duB dT 
sinOsdriu/dT 

1 - sin^ Osu'^ 



ST « 2.8 X 10"^ rad, (23) 



where we have used riD = 1-33, « 10~^ for 
water and we have taken ST = 1 K. In a typical small 
angle setup the sample-to-CCD distance is at least 10 
cm and the pixel size is about 10 /im, resulting in a shift 
of the speckle associated to the direction 9^ of the order 
of 0.28 pixels, sufficient to significantly reduce c/. Thus 
temperature fluctuations of the order of 1 K may lead to 
measurable spurious fluctuations of c/ because of refrac- 
tive effects. 

In most single scattering wide-angle apparatuses the 
sample cell is cylindrical and both the incoming beam 
and the scattered light cross all optical interfaces at nor- 
mal incidence, so that refractive effects are avoided. Nev- 
ertheless, any change of the refractive index due to tem- 
perature variations would still result in a change of the 
speckle pattern. This is because each speckle is associ- 
ated to a well defined value of the scattering vector, q^. 
If njj changes, the scattering angle 9^ corresponding to 

Qs is modified, since 9s = 2arcsin^^^^, where Aq is 
the laser in-vacuo wavelength. Accordingly, the speckle 
pattern is contracted or dilated around the q — direc- 
tion. This is the same effect as the radial shift of the 
speckle pattern when changing the wavelength of the in- 
cident radiation, which was studied in the Seventies |53| | . 
When only a limited portion of the speckle pattern cor- 
responding to a small solid angle is imaged on the CCD 
detector (as it is the case, e.g., for the single scattering 
experiments at = 45 deg reported in this paper), such 
a global contraction or dilation results, locally, in a rigid 
shift of the speckles. The change S9s of 9s in response to 
a temperature fluctuation 6T is 



-2Xoqsdnu/dT 



■.5T'. 



d9s dnxi 
-6.2 X 10"^ rad. 



(24) 



where we have used no = 1.33, Aq = 0.535 /im, qs = 
12.44 /xm-i {9s = 45 deg), ^ « 10"^ R-^ and ST = 1 
K. Taking the sample-to-CCD distance to be 10 cm, we 
find that a fluctuation ST = 1 K would result in a shift 
of the speckle pattern of 6.2 /zm, comparable to the pixel 
size. This would lead to a catastrophic drop of cj, since 
the speckle size is typically of the order of the pixel size. 
Indeed, even a fluctuation ten time smaller, ST — 0.1 
K, would have a measurable impact on c/. Similar ar- 
guments apply also to multiple scattering experiments. 
Note that DWS experiments are more sensitive to vari- 
ations of no than single scattering measurements are, 
because minute changes of q at each scattering event add 
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FIG. 15: Second correlation measured for a perfectly frozen 
speckle pattern. The curves are labelled by r, in sec. Note 
the spurious peaks at At = r, due to the CCD dark noise, as 
discussed in the text. 



up along the photon path, eventually resulting in a sig- 
nificant change of the phase of scattered photons. 

Various strategies are possible to avoid the artifacts 
discussed above, or at least to mitigate their effects on 
TRC data. The impact of beam pointing instability can 
be minimized by reducing as much as possible the light 
path between the laser and the sample, or by delivering 
the beam via a fiber optics. In the latter case, beam 
pointing instability results in fluctuations of the laser-to- 
fiber coupling efficiency and thus of the incident intensity, 
/in- Because c/ is normalized with respect to the instan- 
taneous pixel-averaged intensity (see Eq.QJ, fluctuations 
of /in have little if any effect. Temperature should be 
controlled at least to within 0.1 K and the sample or 
sample holder temperature should be monitored, so that 
any suspect feature in c/ could be compared to the tem- 
perature record. A PIV analysis similar to that presented 
in fig.^lis a useful test to check whether a spurious rigid 
shift of the speckle pattern is at the origin of large drops 
of C/. 

The input from PIV measurements may also be used 
to correct for the effect of a rigid shift Ar of the speckle 
pattern. This could be achieved by constructing a cor- 
rected speckle image, shifting back by — Ar the second 
image of the pair used to compute c/. The corrected 
image would then be used to calculate c/. Standard 
image-processing techniques [H^ can be used to obtain 
sub-pixel shifts. Alternatively, one may exploit the fact 
that in the presence of both a rigid shift Ar and a gen- 
uine evolution of the speckle pattern, the measured de- 
gree of correlation, c/.moas, factorizes as c/^„ioas(^, t) = 
t) [(72.sD(Ar) — 1] (a proof is given in the Appendix 
of ref. [l3|)- Here g2,sp(Ar) — 1 is the normalized spatial 
autocorrelation function of the speckle image and c/(t, t) 
is the degree of correlation that would be measured in the 
absence of any shift. We are currently exploring both ap- 
proaches. 

We conclude this section by reporting a spurious fea- 
ture due to the electronic noise ep, which may affect the 



second correlation Cc, (At, r) measured for time lags r 
much smaller than the typical relaxation time. On those 
short time scales, the speckle pattern is essentially frozen; 
therefore, we shal use as an example the series of im- 
ages of a perfectly frozen speckle pattern obtained as 
described at the end of subsec. IIV CI We process the 
data as usually and calculate the variance of c/. For all 
lags r > the same value is obtained, which is taken as 
an estimate of . This has the advantage of excluding 
other artifacts, such as a possible rigid shift of the speckle 
pattern. As shown in fig. ^1 the second correlation ex- 
hibits spikes at At — r, emerging form a baseline close 
to zero. This is a spurious effect due to the fact that 
the electronic noise Cp is delta-correlated in time. There- 
fore, tp contributes an extra term to c/(t, T)ci{t -\- At, r) 
whenever two out of the four intensity terms involved in 
this expression are measured at the same time, as for 
At = r. A rigorous calculation can be found in ref. [35l |. 
showing that the height of the spurious peaks is 0.5, in 
good agreement with fig. ^1 The spikes of the second 
correlation shown in a figure of a (withdrawn) preprint 
authored by some of us |55J are most likely due to this 
effect. 



VII. CONCLUSIONS 

We have shown that TRC allows heterogeneous dy- 
namics to be unambiguously discriminated from homo- 
geneous dynamics. In order to quantify dynamical het- 
erogeneity, three statistical objects are particularly in- 
sightful: the variance of c/, which corresponds to the dy- 
namical susceptibility X4, the PDF of the degree of corre- 
lation, and its time autocorrelation function -the second 
correlation. Statistical noise due to the finite number of 
pixels over which c/ is averaged can contribute signif- 
icantly to the fluctuations of the degree of correlation. 
Thanks to the correction scheme described in this paper, 
the variance, PDF, and time autocorrelation of c/ can be 
corrected for this contribution, thus making quantitative 
measurements of heterogeneous dynamics accessible to 
scattering techniques. These corrections are particularly 
important when the intrinsic fluctuations are comparable 
to or even smaller than the noise. This may occur be- 
cause dynamic heterogeneity is mild, e.g . in moderately 
concentrated suspensions of colloids [5y| , or because the 
number of available pixels is reduced, e.g. in a small angle 
light scattering or XPCS setup, where the lowest acces- 
sible q vectors correspond to small rings centered around 
the transmitted beam and containing a limited number 
of pixels |20| . Corrections are also important when com- 
paring TRC data obtained from different apparatuses, 
for which both the speckle size and the number of pix- 
els may differ, or when analyzing small angle data, since 
the number of pixels in the rings associated to different 
q vectors varies with q. 

In addition to providing a quantitative description of 
dynamical fluctuations for systems in a stationary or 
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quasi-stationary state, TRC is a useful tool for study- 
ing rapidly evolving dynamics, e.g. during gelation |35| . 
as well as the response to an instantaneous perturbation, 
e.g. an applied shear [H^l- Because in both cases the dy- 
namics may evolve on time scales comparable to or even 
shorter than the relaxation time of g2, a representation 
of the time evolution of c/ is more appropriate and in- 
sightful than that of the two-time correlation function. 
Limitations to the applicability of TRC are mainly due 
to the CCD acquisition rate, which typically does not 
exceed a few tens or hundreds of Hz. An additional ex- 
perimental constraint is the need to store the acquired 
images on the hard disk of a PC and to process them at 
the end of the experiment: data sets up to several Gb are 
not infrequent. 

Other techniques have been proposed in the last few 
years to study dynamical heterogeneity. New light scat- 
tering methods include Speckle Visibility Spectroscopy 
(SVS) |58j, |59i] and the measurement of higher order in- 
tensity correlation functions 0j In a SVS exper- 
iment, one measures ci(t,T = 0), that is the instanta- 
neous contrast — or visibility — of the speckle pattern, 
< Ip{t) >p / < Ip{t) >p — 1. The contrast depends on 
the evolution of the speckle pattern during the exposure 
(integration) time of the CCD. A significant evolution on 
this time scale leads to a blurred speckle pattern image, 
and thus to a reduced contrast. Fluctuations of the con- 
trast can therefore be related to dynamical fluctuations 
on the time scale of the exposure time. Because the latter 
is typically much shorter than the time between succes- 
sive CCD acquisitions, SVS and TRC provide comple- 
mentary information, on fast and slow dynamics, respec- 
tively. Note that SVS data can be obtained on the fly, 
with no need to store images, since only one image at a 
time has to be processed. Higher order intensity corre- 
lation functions 42, J^] , calculated by a dedicated hard- 
ware, allow one to discriminate between homogeneous 
and intermittent dynamics. The time scales that are 
probed and the required measuring time are similar to 
those in a traditional light scattering experiment: dy- 
namics as fast as a fraction of //sec can be measured, but 
the largest available delay is limited to a few tens of sec- 
onds and the experiment duration has to be at least 1000 
times longer than the largest relaxation time of the sys- 
tem. These constraints limit the applicability to glassy 
soft matter. 

Video and confocal scanning microscopy have been 
used to study dynamical heterogeneity in concentrated 
colloidal suspensions '63|. Microscopy and TRC are 
complementary techniques: while the former provides 
unsurpassed details on the motion at a single particle 
level, scattering data typically benefit from better statis- 
tics. Additionally, scattering experiments usually require 
smaller particles than microscopy, a plus when dealing 
with very slow dynamics and when sedimentation effects 
should be minimized. Other techniques that can probe 
heterogeneous behavior in soft glasses are dielectric [gJ 
and rheological |62j | measurements whose sensitivity is 



pushed to the limits set by thermal fluctuations. The 
outcome of these experiments shares intriguing similari- 
ties with TRC data, such as the non-Gaussian distribu- 
tion of voltage fluctuations in dielectric measurements on 
suspensions of Laponite 61J. Whether these similarities 
are coincidental or stem from a common physical origin 
remains an open question. 

In the future, a deeper understanding of both spon- 
taneous dynamical fluctuations and the response to ex- 
ternal perturbations in soft glassy systems will likely re- 
quire the combination of different techniques, possibly 
applied simultaneously on the same sample. We believe 
that TRC can play an important role in this endeavor, 
thanks to its ability to detect instantaneous variations in 
the dynamics. 



Acknowledgments 

We thank P. Ballesta, E. Pitard, L. Berthier, P. 
Holdsworth, J. P. Garrahan for many useful discussions, 
and D. Popovic for bringing to our attention the works on 
the second spectrum. This work is supported in part by 
the European MCRTN "Arrested matter" (MRTN-CT- 
2003-504712) and the NoE "Softcomp", and by ONES 
(grant no. 03/CNES/4800000123), CNRS (PICS no. 
2410), Minstere de la Recherche (ACI JC2076) and the 
Swiss National Science Foundation. L.C. is a junior mem- 
ber of the Institut Universitaire de France, whose support 
is gratefully acknowledged. 



APPENDIX A: CALCULATION OF THE 
VARIANCE AND COVARIANCE TERMS IN 
EQUATION © 

In this section we provide additional details on the two 
approaches to the calculation of the variance of the mea- 
surement noise mentioned in sec. IIIII We first show that 
the variance and covariance terms in Eq. ® may be ex- 
pressed in a way that is independent of the homogeneous 
vs. heterogeneous nature of the dynamics, and then ex- 
plicitly demonstrate the N~-^ scaling of that was used 
in the correction scheme presented in this paper. 

The quantities of interest are the variance of the pixel- 
averaged intensity, aj, the covariance between / and J, 
the variance of the un-normalized intensity cor- 
relation function, aQ^{T), and the covariance between 
G2 and /, (tc2j{t). We recall that I{t) =< Ip{t) >p, 
J{t) =< Ip{t + r) >p and drop in the following the ex- 
plicit dependence on t in the notation of all time- varying 
variables. We assume the dynamics to be stationary and 
the experiment duration much longer than the average 
relaxation time of c/ . We start by noting that cr| quanti- 
fies the fluctuations of the instantaneous pixel-averaged 
intensity. Therefore, ctj is independent of the nature of 
the dynamics: indeed, any dynamical process that fully 
renews the speckle pattern will allow all possible values of 
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TABLE L Coefficients needed to calculate the contribution 
of the measurement noise to iji^j{t), ctq (t), and ac2,iiT) . 





r = 


r ^ 00 

















O"G2{0),/ 





/ to be sampled over time and thus will lead to the same 
aj. In contrast, ai^j{T), aQ^{T), and ctgs, /(■'') depend on 
the way the correlation between distinct speckle patterns 
fluctuates with time and therefore contain contributions 
due both to the measurement noise and to the intrin- 
sic fluctuations of the dynamics. However, note that the 
T — > and T 00 limits of these quantities are insensi- 
tive to the dynamics, because they involve either instan- 
taneous properties of the speckle pattern (for t — 0), or 
quantities related to pairs of speckle images totally un- 
correlated (for r —^ 00). This observation, together with 
the linear dependence of uij(t), (T^^(r), and ug2,i{t) 
on cj in the absence of dynamical heterogeneity shown 
in fig. 121 suggests the following form for the contribu- 
tion of the the measurement noise to the covariance and 
variance terms: 



o'/,j(t) = 0-7,7(00) -I- 



c/(t) 



<^G2,/(t) = crG2,/(oo) 



[o-/,j(0) - cr/,j(oo)] 



+ ^[^G.j(0)-aG.j(oo)]. (Al) 

All coefficients in the r.h.s. of the above equations are 
not affected by dynamical heterogeneity. Table I sum- 
marizes their values in terms of quantities that can be 
directly obtained from the speckle images, regardless of 
the nature of the dynamics. 

As mentioned in sec. IIVI estimates of cr^ obtained by 
using Eq. (|A1|) and Table I are typically affected by a 
significant error. In this paper we have proposed a more 
robust correction method, based on the scaling of 
the measurement noise variance, which we demonstrate 
here. For the sake of simplicity, we assume in the follow- 
ing that there is no correlation between the instantaneous 
value of the intensity at distinct pixels, i.e. Iplq — Ip Iq 
for p ^ q. Physically, this corresponds to a speckle size 
much smaller than the pixel size; this requirement con- 
siderably simplifies the calculations, but it can be relaxed 
without changing the scaling, as we shall show at 
the end of this Appendix. 

For experiments whose duration is much longer than 
the relaxation time of c7, the intensity at any given pixel 
fully fluctuates many times and its probability distribu- 
tion over time is the same as the instantaneous PDF of Ip 
calculated over all N pixels; therefore averages over time 



and over pixels can be swapped. We take advantage of 
this property and of the statistical independence between 



Ip and Iq to write 



I N _ 2 -2 



(A2) 



This expression may be further simplifled by noting that 
moments of the intensity at any pixel are equal to those 
of, let's say, pixel 1: Ip = I^ and Ip = Ii = I. Hence 



N{N-1) 2 y2 _ 1 -2 



(A3) 



Similarly, one obtains 



^2 X! 



IJ 



N TV ' 



(A4) 



P,Q 

^Ifj, IG2T, 



(A5) 



and 



N 



t2 J2 



N 



(A6) 



Equations (|A3(I - (|A6(I show that all the terms in the l.h.s. 
of the expression of a^, Eq. ij^J, are indeed proportional 

to N-\ 

As a final remark, we show that the form of the equa- 
tions derived above is not changed by a short-ranged cor- 
relation between the intensity of distinct pixels, such as 
that typically observed in the CCD speckle images. We 
consider, as an example, the calculation of erf; similar ar- 
guments apply also to the other variance and covariance 
terms. Because the spatial correlation is short-ranged, 
the intensity of pixel p will be correlated only to that 
of a small number, M, of nearby pixels, whereas for all 
other pixels Iplq « Ip Iq- We indicate the set of nearby 
pixels by Qp and split the double sum over distinct pixels 
in Eq. (jA2|) as follows: 



p¥=(i 



P,qeQ 
1 

N 



]Sf2 -P-1 

p,qiQp 



q&Ql 
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By substituting Eq. HA7|I in Eq. (|A2|) . one finds that aj 
still scales as N^^, although with a different prefactor: 



2 1 72 1 



qeQi 



M-2 

TV ■ 



(AS) 



APPENDIX B: PDF OF cj 

We show here that in the limit of large A'' and for 
homogeneous dynamics the PDF of c/ is Gaussian, as 
found experimentally (see fig. El . Additionally, we will 

demonstrate the relationship G2 — I (c7+ 1). Note that 
this relationship is not trivial, since in general (c7+ 1) = 



I. J 



differs from 



<ipJp>p 
77 



G2 



We start by noting that G2 , /, and J are obtained from 
an average over a large number N of pixels. Because of 
the central limit theorem, their PDF is Gaussian, with 
a standard deviation much smaller than the mean. Ac- 



cordingly, at any time G2(i) = 6*2 [1 + eciWii with eq^ a 
Gaussian distributed random variable with mean Tc^ — 
and variance ~ N^^. One can write similar expressions 
for I{t) and J{t), so that to leading order 



G2 
IJ 



G2{l + eG2 



/(l + e/)J(l + ej) 



^(1 
-2(1 



0, 



(Bl) 



where we have used J — I and have introduced ^ — 
£G2 ~ ~ ^ is the sum of three (partially correlated) 
Gaussian random variables and therefore is itself Gaus- 
sian distributed 47], with mean — eJ ~ ej = 0. It 
follows that the PDF of 



(B2) 



—2 , 



is Gaussian, with mean c/ = ^ — 1. Thus, G2 = I (ci - 
!)• 
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